Quantum effects in a weakly-frustrated S=l/2 two-dimensional Heisenberg 
antiferromagnet in an applied magnetic field 
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We have studied the two-dimensional S=l/2 square-lattice antiferromagnet Cu(pz)2 (0104)2 using 
neutron inelastic scattering and series expansion calculations. We show that the presence of anti- 
ferromagnetic next-nearest neighbor interactions enhances quantum fluctuations associated with 
resonating valence bonds. Intermediate magnetic fields lead to a selective tuning of resonating va- 
lence bonds and a spectacular inversion of the zone-boundary dispersion, providing novel insight 
into 2D antiferromagnetism in the quantum limit. 
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Some of the best examples of macroscopic quantum 
magnets can be found in low-dimensional antiferromag- 
nets that preserve strong quantum fluctuations to very 
low temperatures. Cases in point are antifcrromagnetic 
(AF) chains that are quantum critical at zero tempera- 
ture and feature deconfined spinon excitations for S=l/2 
but a spin liquid with gapped triplet excitations for 
S—1 In two dimensions, quantum effects are gener- 
ally reduced and the ground state of S=l/2 square lat- 
tice Heisenberg antiferromagnets is long-range ordered, 
albeit with a reduced ordered moment and spin-waves 
whose energy is renormalized by quantum fluctuations 
[^, Ql . Importantly, however, short range quantum cor- 
relations are preserved in the ground state and lead to 
a quantum-induced dispersion at the AF zone bound- 
ary and to a weak continuum of states at high energies 

I,i,i0,ii[i3,[iil,[il. 

Competing interactions generally increase quantum 
fluctuations, particularly for materials close to a quan- 
tum critical point. For the 2D S=l/2 antiferromagnet, 
one source of competition can be next-nearest neighbor 
interactions that destabilize a simple nearest-neighbor 
antiparallel alignment and thus enhance fluctuations 
[ll,ll3|. Such a scenario has been extensively investi- 
gated using the Ji — J2 model, where Ji and J2 give the 
stren gth of the nearest and next-nearest neighbor interac- 
tions [13| . The Ji — J2 model features two different types 
of order for J2/J1 < 0.38 and J2/J1 > 0.6, and proba- 
bly a spin-liquid phase in between. The behavior of the 
2D S=l/2 antiferromagnet in magnetic fields was studied 
theoretically for the case of nearest nei ghb or interactions 
[3, El, m as well as for Ji - J2 model 0, Ell . However, 
little is known experimentally about the magnetic field 



behavior of the 2D S=l/2 antiferromagnet, particularly 
in the presence of next-nearest neighbor interactions. 

We have studied a weakly-frustrated 8=1/2 AF square- 
lattice where next-nearest neighbor interactions are a 
small perturbation that allows us to study the change 
of the magnetic properties upon introduction of com- 
peting interactions. Our inelastic neutron scattering ex- 
periments were performed using the 2D square lattice 
Heisenberg antiferromagnet Cu(pz)2 (0104)2. Deuter- 
ated Ou(pz)2(C104)2 crystallizes in the monoclinic C2/c 
space group [20|. Each Ou^"*" ion in Cu(pz)2 (0104)2 car- 
ries S=l/2 and is surrounded by two pairs of identical cis 
pyrazine molecules, creating a two dimensional square 
array of copper atoms linked by pyrazine molecules in 
the be plane. The two fold rotation axis (0, y, 1/2) and 
the mirror plane parallel to the ac plane ensures that all 
nearest neighbor interactions between Cu^^ ions on the 
square lattice are identical. From the temperature de- 
pendence of the magnetic susceptibility, it was concluded 
that Cu(pz)2 (0104)2 represents a S=l/2 2D square lat- 
tice antiferromagnet with nearest-neighbor exchange J = 
1.53(8) meV and a saturation field of /io-ffsat ~ 45 T 



22l |. Due to interplane interactions, Ou(pz)2 (0104)2 

adopts long-range AF order below Tn = 4.2 K and the 

ratio between the interlayer to intr alay er exchange has 

4 



201 2 



been estimated as J'/ J — 6.8 • 10" 

The neutron experiments were performed using the 
cold-neutron triple-axis spectrometers PANDA at FRM2 
at Oarching and IN14 at ILL at Grenoble, and an ar- 
ray of deuterated Ou(pz)2 (0104)2 single crystals with a 
total mass up to 2g co-aligned with a mosaic of 0.5°. 
The sample was aligned with its reciprocal (0, k, 1) 
plane in the horizontal scattering plane. A cryomag- 
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FIG. 1: Zone boundary spin dispersion at zero field and 
H ~ J. Color plot of the normalized scattering intensity 
J(Q,tj) at T = 80 mK, showing the dispersion from (vr, 0) 
to (7r/2, 7r/2) at zero field and /xq// = 14.9 T under otherwise 
identical conditions. Both panels present smoothed data ob- 
tained by performing six constant Q-scans from (7r/2, 7r/2) to 
(tt, tt) with a 0.05 meV energy step. 



net allowed the application of vertical magnetic fields 
up to 14. 9T. The magnetic field was thus nearly per- 
pendicular to the square-lattice plane. To probe the 
ground state as a function of magnetic field, we used 
a dilution refrigerator, reaching temperatures of the or- 
der of 50-lOOmK in each of the two experiments. The 
measurements were performed using a fixed final energy 
Ei = 4.66 meV and E{ = 2.98 meV for the PANDA 
and IN14 experiments, respectively, obtained via the 
(002) Bragg reflection from a pyrolithic graphite (PG) 
monochromator, a focused analyzer and a cooled Be fil- 
ter before the analyzer. The chemical unit cell contains 
four Cu^"*" atoms, thus nearest-neighbor AF order does 
not break translational invariance, and Q — (0, 1,0) and 
Q = (0,0, 1) correspond to the AF point (7r,7r). Conse- 
quently, Q = (0,TO ± 1/2, n± 1/2), where m and n are 
integers, correspond to (tt, 0), while Q = (0, m ± 1/2, n) 
corresponds to (7r/2, 7r/2). 

A color plot of the normalized neutron scattering spec- 
tra at the zone boundary is shown in Fig. [T] for zero ap- 
plied field and hqH = 14.9 T. At zero field, the onset 
of scattering at reciprocal wave-vector (tt, 0) is reduced 
by 11.5(7)% in energy compared to (7r/2,7r/2), at odds 
with spin-wave theory. This is larger than expected from 
Quantum Monte Carlo simulations and series expansion 
calculations for the Heisenberg square-lattice antiferro- 
magnet with nearest-neighbor interactions, which predict 
a zone-boundary dispersion of 8% to 10% [3,0]. Our ob- 
servation is opposite to what has been observed in the 
high-Tc material La2Cu04 where the energy at (tt, 0) is 
higher than at (7r/2,7r/2) [23|. The latter has been at- 
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FIG. 2: Series expansion calculations of the zone-boundary 
dispersion. Theoretical magnon dispersion for diff'erent values 
of the next-nearest neighbor exchange interaction J2 in zero 
magnetic field (left panel) and in a fioH = 14.9 T magnetic 
field (right panel) . The value of the nearest neighbor interac- 
tion J has been normalized to fit the experimental data. The 
normalization constants are shown in the legend. 



tributed to ring-exchanges arising from fmite-U/t 24 1. 

We developed series expansion calculations for the 
Ji — J2 model in a magnetic field 25, 2^, 2^. Since the 
magnetic field causes the spins to become non-collinear, 
this requires the technically difficult multi-block method 
to be used [S^l . Our series expansion calculations for the 
Ji - J2 model [H, [ii, [23 show that the addition of J2 
leads to an enhancement of the zone boundary disper- 
sion. Fig. [5] shows the calculated single- magnon energies 
for different relative strengths J2/J1 between nearest to 
next-nearest neighbor interactions. For Cu(pz)2 (0104)2 
we thus attribute the increased zone-boundary disper- 
sion to a small AF next-nearest neighbor interaction of 
the order of J2 - 0.02 - 0.05 J. 

According to spin wave theory, the energy at {tt/2, 7r/2) 
is equal to -E(7r/2,7r/2) — 2(Ji — J2) at all fields. Next- 
nearest neighbor interactions thus lead to a smaller 
zone-boundary energy, and thus to a smaller effective 
nearest-neighbor exchange J. Using J = Ji — J2 and 
J — 1.53(8) meV determined from susceptibility mea- 
surements [2^ [13], the field-induced change of £^(7r/2,7r/2) 
is obtained from the field dependence of the renormaliza- 
tion factor Zc = -E(7r/2,7r/2)/2J^, yielding Zc = 1.19(2) at 
zero field in excellent agreement the predicted value of 
Z, - 1.18 0. 

Fig. [3^ shows that the spin excitation at (tt, 0) is 
considerably broader than experimental resolution, and 
that neutron scattering extends to about E = 3.9 meV. 
This represents clear experimental evidence of a mag- 
netic continuum at (tt, 0) of the square-lattice antiferro- 
magnet that is not expected from spin-wave theory. The 
continuum is clearly stronger than in copper deutero- 
formate tetradeuterate (CFDT) with a zone-boundary 
dispersion of 7% where next-nearest neighbor exchange 
are absent suggesting that next-nearest neighbor ex- 
change enhances continuum excitations in the 2D S=l/2 
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FIG. 3: Field dependence of zone-boundary excitations. En- 
ergy scans at (tt, 0) and (7r/2, 7r/2) are shown in (a) and (b) re- 
spectively. The data was measured at T = 80 mK. The black 
solid line in (b) represents linear spin- wave excitations convo- 
luted with the resolution function. Apart from the continuum 
region, all the experimentally measured peaks are resolution 
limited. The field dependent onset of magnetic scattering at 
(tt, 0) and (7r/2, 7r/2) as a function of energy is shown in (c) 
by red triangles and black circles respectively. The dashed 
red and black lines represent the estimate from the series ex- 
pansion calculations with Ji — 1.54 meV and J2 = 0.02Ji. 
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FIG. 4: Spin dispersion for H ~ J. (a) The spin wave dis- 
persion in Cu(pz)2 (0104)2 measured at fioH = 12 T and 
T = 80 mK. The red line represents the spin- wave dispersion 
with Ji = 1.54 meV and J2 — 0.02Ji. (b,c) Energy scans at 
(O.Gtt, O.Gtt), (O.Ttt, O.Ttt) provide evidence of a mode with a 
field-induced gap. (d) Constant-energy scattering at 1 meV 
energy transfer providing evidence for the Goldstone mode. 
The solid line in (b)-(d) corresponds to a convolution of a 
Gaussian with the resolution function, demonstrating that 
the excitations are resolution limited. 



square-lattice antiferromagnet. The extended continuum 
in Cu(pz)2 (0104)2 is also in contrast to well-defined exci- 
tations observed in La20u04 near (tt, 0), suggesting that 
ring exchange has the opposite effect on the continuum 
to next-nearest neighbor interactions. Further, Fig. [1] 
also provides evidence of an energy gap around {tt/2, n/2) 
that separates a main mode and a much weaker contin- 
uum above 4.2 meV, as predicted by Ho et al. 

Fig. [3]; shows that magnetic fields strongly affect the 
quantum fluctuations at the zone boundary: The energy 
of the excitation at (7r/2, 7r/2) decreases much faster with 
field than that at (tt, 0). The zone-boundary dispersion 
at /io-ff = 14.9 T is inverted from that at zero field. 
Our series expansion calculations show that, with the 
application of a magnetic field, the inversion of the zone 
boundary dispersion only occurs for sufficiently small 
J2/J1. A five percent J2/J1 no longer shows the rever- 
sal of the dispersion that is seen in experiments. This 
implies that even a relatively small next-nearest neigh- 
bor interactions are effective in enhancing the contin- 
uum of excitations, consistent with our estimate that 
J2/J1 = 0.02. The quantum renormalization factor Zc 
decreases rapidly from Zc = 1.19(2) at zero field and ap- 
proaches Zc — 0.99(2) at hqH = 14.9 T, also consistent 
with our calculations. 

It is known that the wave-vector dependence of the 
magnetic excitations at the zone boundary is a result of 
resonating valence bond quantum fluctuations between 
nearest-neighbor spins that reduce the energy of resonant 



excitations near (tt, 0) below what is expected from renor- 
malized spin- wave theory 0, |^ . The field- induced rever- 
sal of the zone-boundary dispersion reveals that mag- 
netic fields of the order of 77 ~ J strongly couple to 
these local quantum fluctuations. This dispersion, that 
is not expected from spin-wave theory, demonstrates the 
presence of local quantum fluctuations in the 2D S=l/2 
square-lattice antiferromagnet even for fields of the order 
of -ff ~ J. The quantum origin of the dispersion rever- 
sal is also confirmed by our series expansion calculations. 
The observed dispersion suggests that the energy of the 
resonance at (7r,0) is raised by 4.5(7)% above that of 
renormalized spin-wave theory, providing direct evidence 
of field-tuned resonating valence bond fluctuations. 

The excitation spectrum at fioH = 12 T features two 
well defined magnetic modes of excitations (Fig. |3|). The 
spectrum consists of a Goldstone mode, that indicates 
unbroken rotational spin symmetry in the plane perpen- 
dicular to the magnetic field, and a gapped mode. The 
field dependence of the mass of gapped mode at the AF 
zone center, shown in Fig. [SI is linear, which is consistent 
with the theoretical arguments discussed in [2^. How- 
ever the finite gap at zero field, obtained by the linear fit 
to the experimental data, implies the existence of a small 
anisotropy in the system and requires additional inves- 
tigations. We also analyzed the spectral weight of the 
gapped excitation at (tt, tt). Our calculations completely 
describe the observed field dependence of the peak inten- 
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FIG. 5: Field dependence of the zone-center excitation. En- 
ergy scans performed at (7r,7r) at ^oH = 2 T, fioH = 6 T, 
HqH = 10 T (a) and the gap energy at the AF point plot- 
ted as the function of field (b). The data was measured at 
T — 80 mK. The solid lines in (a) are the fits of a Gaussian 
function convoluted with the resolution function. The filled 
circles in (b) represent the experimental data. The solid line 
in (b) is the linear fit of the gap energy. Black circles in in- 
set (c) show the measured intensities as function of magnetic 
field and the curve is the scattering intensity calculated using 
linear spin-wave theory. 



sities of the gapped mode at the AF zone center (Fig.[5j;). 

Fig. m shows the observed dispersion from (tt, tt) to 
(7r/2, 7r/2) compared to hnear spin-wave theory. The hn- 
ear spin-wave theory calculations in a magnetic field was 
performed for the following Hamiltonian: 



n = .h }_^s, 

<i,j> 



J2 S.-Sfe-A^oi^E^?' 

<i,k> i 



where represents nearest neighbor pairs of spins, 

(j, k) next nearest neighbor pairs and is a spin com- 
ponent perpendicular to the be plane. The renormal- 
ized spin wave theory with Joff — 1.54(l)meV and Zc — 
1.03(1) for ^qH = 12 T describes quahtatively the ob- 
served dispersion, as shown in Fig. |4ja). However, for 
wave- vectors near (37r/4, 37r/4), spin- wave theory clearly 
predicts a larger spin-wave energy than experimentally 
observed. This difference may reflect the importance of 

Possibly, better 
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magnon-magnon interactions 
agreement could be obtained by including higher-order 
terms in the spin-wave calculation. 

In summary, our study shows that the zone-boundary 
dispersion and the zone-boundary continuum of the 2D 
S=l/2 square-lattice antiferromagnet are enhanced by 
next-nearest neighbor interactions, while the spin-wave 
energies merely experience a small renormalization. Mag- 
netic fields of the order of H ^ J lead to a qualitative 
change of the quantum fluctuations that suppress the 
continuum of excitations and renormalize the spin-wave 
velocity, but without suppressing the zone-boundary dis- 
persion that arises from non-trivial quantum fluctua- 
tions. In fact, in contrast to spin- wave theory, we find 



that the zone-boundary dispersion is inverted compared 
to zero field, providing direct evidence of a field-induced 
change of resonating valence bond fluctuations. 
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After submission of our manuscript, we became aware 
of A. Luscher and A. Laeuchli [29!| who also find theo- 
retically that the zone boundary dispersion inverts as a 
function of field. 
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